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Abstract 

We use the tools of hybrid intranuclear-cascade/nuclear-de-excitation models to evaluate the 
sensitivity of several physical observables to the inclusion of a multifragmentation stage in the 
de-excitation chain and assess the need for a multifragmentation model in the quantitative de- 
scription of p+^^Fe and p+^^^Xe reactions at 1-GeV incident energy. We seek clear signatures 
of multifragmentation by comparing different state-of-the-art de-excitation models coupled with 
intranuclear-cascade models and by focusing on discriminating observables such as correlations and 
fragment longitudinal- velocity distributions. None of the considered observables can be unambigu- 
ously interpreted as a multifragmentation footprint. The experimental data are best described 
as originating from sequential binary decays. However, no de-excitation model can reproduce the 
experimental longitudinal-velocity distributions from 1-GeV p+^'^^X.e. 
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I. INTRODUCTION 



Multifragmentation is generally considered as the quasi-simultaneous break-up of highly- 
excited nuclear matter into clusters and unbound nucleons. Interest towards this phe- 
nomenon was first triggered by anomalously large production cross sections of intermediate- 
mass fragments (abbreviated as IMFs and defined as 3 < Z < 10 for the purpose of this 
paper) from intermediate-energy heavy-ion collisions (see Refill for a collection of recent 
reviews). The earliest theoretical explanations suggested to interpret the typical power-law 
distribution of fragment masses as a signature of liquid- vapour equilibrium of nuclear matter 
near the critical temperature. In this scheme, nuclear clusters are simultaneously formed, 
with large multiplicities, as condensation droplets of a vapour of nucleons. However, as the 
initial enthusiasm over liquid-vapour multifragmentation faded and other candidate models 
(e.g. statistical multifragmentation, spinodal instability and even sequential binary decays) 
were put forward to explain the data, it was quickly realised that the power-law signature 
was by no means unique of liquid-vapour multifragmentation. Remarkably, even simple 
percolation models are able to reproduce most of the features of the observed IMF distribu- 
tions. Therefore, other observables must be sought if one wishes to discriminate among the 
proposed IMF production mechanisms, which are anyway not necessarily mutually exclusive. 

One of the main difficulties of multifragmentation studies based on heavy-ion reactions 
is that there is considerable theoretical uncertainty on the reaction dynamics and on the 
importance of collective effects such as deformation or compression. In nucleon-induced 
reactions, on the other hand, it is difficult to imagine that the collective state of the system 
can be strongly perturbed. Since it had been known for a long time that IMFs could also 
be produced in nucleon-induced reactions, multifragmentation studies were also performed 
on these better-understood systems, although the excitation energies that can be reached 
by this method are typically lower. 

Today, the importance of multifragmentation in nucleon-induced reactions is the subject 
of a long-standing discussion. While it is generally accepted that multifragmentation will 
eventually set in at high projectile energy, due to the increasing energy transfer from the 
projectile to the target nucleus, it is not yet clear whether and to what extent multifrag- 
mentation needs to be postulated for a reliable quantitative description of reactions around 
1 GeV, a region which is most interesting for technical applications such as Accelerator- 
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Driven Systems (ADS) 3, radioprotection in space [3|| and shielding at accelerators 

The recent IAEA-promoted "Benchmark of Spallation Models", which focused on the 60- 
3000-MeV incident-energy range, represents an effort "to assess the prediction capabilities 
of the spallation models used or that could be used in the future in high- energy transport 
codes; to understand the reason for the success or deficiency of the models in the different 
mass and energy regions or for the different exit channels; to reach a consensus, if possible, 
on some of the physics ingredients that should be used in the models" j^. The bench- 
mark saw the participation of seventeen spallation models, all of which were couplings of 
a dynamical reaction model (intranuclear cascade, quantum molecular dynamics. . . ) and a 
statistical-decay model, with the possible presence of an intermediate pre-equilibrium stage. 
Since not all the participating models include a multifragmentation stage, it is in princi- 
ple possible to study the benchmark results and estimate the sensitivity of the benchmark 
endpoints (isotopic production cross sections, excitation curves, neutron-multiplicity distri- 
butions, double-differential cross sections for neutrons, light charged particles and pions) 
to the multifragmentation process. In particular, by comparing the predictions of differ- 
ent de-excitation models coupled with a fixed dynamical stage, one can extract precious 
information about the influence of de-excitation alone. 

However, previous studies have already indicated that inclusive observables, such as 
double-differential nucleon spectra or nuclide yields, are rather insensitive to the inclusion 
of a multifragmentation stage in the de-excitation chain Hence, characteristic signatures 
of multifragmentation must be sought among other, more discriminating observables. The 
impact of a multifragmentation stage in the de-excitation chain can in principle be assessed 
by comparing calculation results with experimental data. 

The goal of the present work is to identify possible signatures of multifragmentation 
by studying nucleon-induced reactions with the tools of coupled intranuclear-cascade/de- 
excitation models. We shall focus on the 1-GeV p+^^Fe and ^+^^^Xe reactions, which 



have been the object of recent studies 6l-|l0|. The small mass of the p+^^Fe system leads 



to the production of a limited number of nuclides. Several de-excitation mechanisms can 
contribute to a given nuclide yield, making it more difficult to extract an unambiguous 
multifragmentation signature from a background of de-excitation residues and/or direct 
IMF emissions. On the other hand, the multifragmentation threshold may be more easily 
attained in p+^^Fe, which realises higher excitation energies per nucleon. Thus, the two 



systems studied are complementary. Heavier systems are excluded from the present study, 
in order to avoid the conceptual and technical complications connected to the competition 
between light-fragment emission and fission. 

The paper is organised as follows: in Sec. [Til we give a brief overview of the models 
used for the study. Sec. IIIII examines the sensitivity of the considered observables to the 
choice of the intranuclear-cascade model. Sec. IIVI discusses the available inclusive residue- 
production data and how they are reproduced by the different de-excitation models. Sec. IVl 
presents the model predictions for the SPALADIN correlation data-set for p+^^Fe 6|. Sec. lVIl 
discusses longitudinal-velocity distributions measured in inverse kinematics at GSI, Darm- 
stadt, Germany ^, [lo|. Sec. IVIII discusses time intervals between fragment emissions in 
INCL4.5/GEMINI++. Finally, Sec. Ivml summarises our conclusions. 



II. MODEL OVERVIEW 



All the calculations presented in this paper were performed using a coupled intranuclear- 
cascade/statistical-de-excitation model. We used two intranuclear-cascade model (INCL4.5 
and Isabel) coupled with three different de-excitation models (ABLA07, GEMINI++ and SMM). 
Since the focus of this paper is on de-excitation, we will limit ourselves to directing the reader 
to the relevant publications for details about the physics of the cascade models. 



A. Cascade models 

The Liege Intranuclear Cascade model (INCL) is one of the most refined existing 

tools for the description of nucleon-, pion- and light-ion-induced reactions in the 150-3000- 
MeV incident energy range. The model is currently maintained and developed jointly by 
the University of Liege (Liege, Belgium) and CEA (Saclay, France). It can describe the 
emission of nucleons and pions; light clusters (up to Z = 5, A = 8 with the default program 
options) can also be produced through a dynamical phase-space coalescence algorithm. The 
INCL model is not to be considered as adjustable. It does contain parameters, but they are 
either taken from known phenomenology (such as the matter density radius of the nuclei) 
or have been adjusted once for all (such as the parameters of the Pauli blocking or those 
that determine the coalescence module for the production of the light charged clusters). 
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The predictions of INCL concerning those observables that can be confronted directly to 
experiment, namely the high energy parts of particle spectra, are of rather good quality, as 
it was recently shown The INCL/ABLA07, INCL/GEMINI++ and INCL/SMM combinations 
were also recognised among the best-performing participants of the IAEA "Benchmark of 
Spallation Models" jsl- The present work is based on the INCL version that was used for the 
IAEA benchmark, plus some minor bug fixes; this version is known as INCL4.5. 

An older version of the INCL model, known as INCL4.2, was employed for studying the 
SPALADIN correlation data-set j^. The most important differences between INCL4.2 and 



4 . 5 are reviewed in Ref. 



12 



and include the introduction of the cluster-coalescence algorithm, 



energy- and isospin-dependent potentials for nucleons and pion potentials, as well as an 
improvement of Pauli blocking. More details are given in the reference above. 

The Isabel model 13|, no longer developed, has vastly contributed to the understanding 
of nucleon-nucleus and nucleus-nucleus reactions. It is mainly used in the present work to 
highlight the sensitivity of the studied observables to cascade. 

B. De-excitation models 
1. ABLA07 

The ABLA07 model {l^ is maintained and developed by the CHARMS group at GSI, 
Darmstadt, Germany. The model contains a multifragmentation sub-module, which is trig- 
gered only if the temperature of the compound nucleus exceeds a mass-dependent (as sug- 



gested by Natowitz et al. 



15|) freeze-out threshold: 



7freeze-out(^) = Hiax (5.5, 9.337 ■ exp(-2.82 ■ 10'^ A)) MeV. (1) 

In that case, the system breaks up into fragments whose mass is distributed according to 
an empirical power-law spectrum, and whose momenta carry Goldhaber-type and thermal 
contributions. Coulomb repulsion among multifragmentation products is accounted for in 
a simplified manner. The excitation energies of the resulting fragments are determined by 
assuming thermal equilibrium at the freeze-out temperature. Subsequent de-excitation of 
the multifragmentation products is assumed to be purely binary. If the multifragmentation 
module is not triggered, the initial compound nucleus directly enters the secondary de- 
excitation phase. 



During secondary de-excitation, emission of any stable nucleus up to half the mass of 
the compound nucleus is possible, and it is quantitatively described by the Weisskopf-Ewing 



evaporation formalism 



16|. Above the Businaro-Gallone point, competition with fission is 



treated dynamically and it is based on solutions of the Fokker-Planck equation for collective 
deformation of the nucleus over the fission barrier. ABLA07's fission module is among the 
most sophisticated models available on the market, but it is only of marginal interest for 
the systems studied in this paper. 

Finally, subsequent binary decays are assumed to be independent; in particular. Coulomb 
interactions among particles produced in different decays are neglected. This assumption is 
customary in binary de-excitation models. 

A less-sophisticated version of the ABLA model was considered for the study of the SPAL- 
ADIN correlations ^ and was found unable to reproduce the measured residue-production 
cross sections in p+^^Fe. It did not include multifragmentation, nor evaporation of frag- 
ments heavier than alpha particles. More details about the differences between the two 



versions can be found in Ref. 
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2. GEMINI++ 



The GEMINI++ model, developed by R. J. Charity 17|, represents an effort to describe 
nuclear de-excitation uniquely in terms of binary decays. No simultaneous break-up is al- 
lowed. Multi-fragment events can of course be produced by sequences of binary fragment 
emissions; as in the case of ABLA07, Coulomb interactions among particles emitted in differ- 
ent decays are neglected. Emission of light particles (Z < 3 by default) is described by the 



lauser-Feshbach evaporation formalism 



181 : Moretto's conditional-saddle-point formalism 



19| with Sierk's finite-range barriers |20j is used for complex-fragment emission. For heavy 



2l|. 



systems, the fission width is calculated using a refined Bohr- Wheeler approach 

GEMINI++'s asymmetric-fission module has recently been improved j22| to describe frag- 
ment yields from fusion and spallation reactions with the same parameter set. To this end, 
it was necessary to augment Sierk's barriers by a constant shift of 7 MeV, which can be 
interpreted as the difference in Wigner energy between the mother nucleus and the nascent 
fragments. However, this interpretation is not devoid of complications; see the relevant 
papers for more details. 
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3. SMM 



The Statistical Multifragmentation Model (SMM) 23|,|2J|, presently maintained by A. S. Botv- 
ina, is one of the most successful and widely applied multifragmentation models. Like in 
ABLA07 the first possible decay stage is the simultaneous break-up of the thermalised source 
in a number of hot fragments and particles. Unlike ABLA07's semi-empirical approach, 
SMM always enters this sub-module and samples break-up configurations according to their 
thermodynamical weight in a given freeze-out volume (taken to be three times the normal 
nuclear volume). Thus, at low excitation energy, single-fragment configurations (i.e. com- 
pound nucleus) naturally dominate the multifragmentation phase; the importance of multi- 
fragment configurations smoothly increases and starts dominating the thermodynamical 
weight around 3 AMeV excitation energy. Mass, charge, excitation energy and momentum 
of the emerging hot fragments are sampled respecting conservation laws. Coulomb acceler- 
ation is then accounted for by solving the Hamilton equations for the propagation of the 
fragments in their mutual Coulomb field. 

Secondary de-excitation is then applied to the hot fragments. If they are sufficiently 
light {A < 16), the Fermi break-up model is applied. Otherwise, according to a modified 
Weisskopf-Ewing jfel] scheme, they can evaporate particles up to ^^O. Fission is described 



by the Bohr- Wheeler model 25 1. 

Note that the SMM version used for the present work employs slightly different evaporation 
barriers compared to the IAEA benchmark. Barriers are computed using the standard 
formula 

B = (1.44 MeV ■ fm) ^^^^ 



In the IAEA benchmark, Tq = 1.5 fm was used. In the present work, Tq is determined as 

1 + 6.103-10-3^1^2 , 
^° = '-'^'- 1 + 9.443. 10-3Ziz/"^- 

This difference is marginal as far as IMF cross sections are concerned. 



III. BASIC CASCADE RESULTS 

As a first, basic comparison. Fig. [1] shows distributions of excitation energy and mass loss 
of cascade remnants of 1-GeV p+^^Fe and p+^^^Xe. Already at this stage, it is possible to 
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Figure 1. Comparison of excitation-energy (left panel) and mass-loss (right panel) distributions for 
remnants of the 1-GeV p+^^Fe and p+^^^Xe reactions, as calculated by different cascade models. 
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Table I. Reaction cross section, average remnant excitation energies (total and per nucleon) and 
average remnant mass predicted by INCL4.5, Isabel, INCL4.5 without cluster coalescence, and 
INCL4.2 for 1-GeV proton-induced reactions on ^^Fe. 



observe that INCL4.5 and Isabel are not equivalent. INCL4.5 produces on average hotter 
and lighter remnants than Isabel, although differences in excitation energy are smaller for 
^^^Xe than for ^^Fe, as quantified in Tables [T] and [Tll Note that the reaction cross sections 
predicte bdy the two codes differ by only a few percent; thus, differences in the remnant 
characteristics must trace back to different cascade histories. 

Previous investigations had found that INCL's and Isabel's excitation-energy distribu- 



tions 



'or p+^^Fe were remarkably similar jo]. The claim concerned version 4.2 of the INCL 



code Uj, which is represented in Fig. [U by the blue lines. This state of affairs was evi- 



dently modified by later developments of the INCL code [12]. We stress that the similarity 
between INCL4.2 and Isabel's results can at least partly be explained by the very similar 
physics content of the models. Excitation energies are also sensibly larger in INCL4.5 than 
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Table II. Same as Table H for 1-GeV p+^^^Xe. 



in INCL4.2. However, this difference cannot be simply ascribed to a single cause, but rather 
represents the combined effect of several new physics ingredients, such as energy-dependent 
nucleon potentials and pion potentials. 

One major difference between Isabel and INCL4.5 (and between INCL4.2 and INCL4.5) 
is INCL4.5's ability to dynamically produce light charged composite particles (see Sec. Ill Al) . 
Fig. [1] depicts INCL4.5's predictions with (solid black lines) and without (dotted black lines) 
cluster coalescence. One immediately observes that Isabel's mass distributions are remark- 
ably similar to those predicted by INCL4.5 without coalescence. Clearly INCL4.5's cluster 
emission algorithm reduces the remnant mass, but does not sensibly affect the excitation- 
energy distribution. This is quite well understood: escaping clusters typically extract spec- 
tator nucleons close to the Fermi sea, thereby reducing the mass of the remnant but without 
significantly affecting its excitation energy. This result is also consistent with the docu- 



mented behaviour of an older, less-refined clustering algorithm 



26|. 



Therefore, we conclude this section by observing that, contrary to what was claimed in 



Ref. |6|, the choice of the intranuclear-cascade model does have some importance. In what 
follows, we shall discuss how the differences in remnant distributions are reflected in the 
residue-production cross sections, and how these observables can guide us in selecting the 
cascade model that should be used for the study of more discriminating observables. 



IV. RESIDUE-PRODUCTION CROSS SECTIONS 



Figures 12] and E] show the predicted residue-production cross sections as a function of the 
nuclide charge, compared with the measurements obtained with the SPALADIN apparatus 
jol or at the FRagment Separator (FRS) j7-9|. The two experiments have different accep- 
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Figure 2. Inclusive residue-production cross sections for 1-GeV p+^^Fe, as a function of the nuclide 
charge. Left panel:^ntranuclear cascade simulated by INCL4 . 5. Right panel: Isabel. Experimental 
data from Refs. 
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Figure 3. Same as Fig. [2l for 1-GeV p+^^^Xe. Experimental data from Ref. [8 



tance cuts. The SPALADIN data, which will be discussed in detail in Sec. |Vl present a 
kinematical cut on particles with large longitudinal velocities with respect to the initial ^^Fe 
nucleus, which are mostly nucleons and light particles. The acceptance was estimated to be 
virtually complete for Z > 4. In the FRS data, on the other hand, only a selected number of 
isotopes were measured. For the model curves in Figs. [2] and [3l we have chosen to define the 
residue-production cross sections as the sums of the calculated isotopic residue-production 
cross sections over the nuclides observed in the FRS experiments. No kinematical cut was 
applied. Since the SPALADIN and FRS data sets are largely compatible, as it clearly 
appears from Fig. |2l we do not expect this choice to bias our analysis. 



10 



For a given de-excitation stage, the curves reflect the differences in the cascade output. 
Note that calculated cross sections for residues close to the target (say Z >20 for ^^Fe and 
Z > 45 for ^'^^Xe) are almost independent of the choice of the de-excitation model and are 
dominated by the cascade model; in particular, they are consistently better reproduced by 
INCL4.5 than Isabel. We consider that an accurate prediction of these cross sections is a 
crucial prerequisite that the cascade model must satisfy if it is to be used for the study of 
more exclusive and discriminating observables. Thus, in Sec. |V] and following, we will only 
retain the INCL4 . 5 model for the analysis of fragment correlations and velocity distributions. 

INCL4 . 5's hotter and lighter remnants also lead on average to lighter residues. Note that 
the IMF cross sections in p+^^Fe and the 10 < Z < 30 cross-section plateau in p+^^^Xe are 
better reproduced by the INCL4.5/de-excitation combinations, while Isabel consistently 
underestimates the ^^^Xe cross sections by a factor of about three. Comparison with the 
excitation-energy distributions in Fig. [1] suggests that these cross sections are associated 
with highly-excited remnants. Likewise, Isabel consistently underestimates cross sections 
with 10 < Z < 18 in p+^^Fe. 

The sensitivity of the inclusive residue-production cross sections to cascade can be fur- 



ther illustrated by considering the results obtained with the GEM de-excitation model |27| . 
coupled with INCL4.5 (left panes of Fig. |2] and ED . In a previous study je], the INCL4.2/GEM 
combination was excluded from the study of SPALADIN correlations because it was un- 
able to reproduce the IMF-production cross sections in p+^^Fe. New INCL4.5/GEM calcula- 
tions predict IMF cross sections that are about a factor of 3 higher than the INCL4.2/GEM 
and in acceptable agreement with the experimental data; this is due to INCL4 . 5's different 
excitation-energy and remnant-mass distributions. However, the plateau cross sections in 
^_l_i36^g are underestimated by at least three orders of magnitude by INCL4.5/GEM. Thus, 
we also exclude the GEM de-excitation model from this study. 

If we now focus on a fixed cascade model (e.g. INCL4.5), we can observe that the three 
de-excitation models produce similar charge distributions. In this sense, we confirm that 
residue-production cross sections are rather insensitive to the de-excitation mechanism. 
However, we remark that de-excitation models present free parameters that can be ad- 
justed to help reproduce the residue-production cross sections. The p+^^Fe data-set, in 
particular, is a popular benchmark for spallation models (cascade/de-excitation) due to 
its good accuracy. The GEMINI++ parameters connected with asymmetric fission were ad- 
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mittedly fitted to the p+ Fe and p+ 'X.e residue-production cross sections, among other 



data-sets 



22|. Thus, Figs. [2] and [3] can deceptively lead to underestimate the sensitivity of 



residue-production cross sections to the de-excitation model. 

The sensitivity to de-excitation can be further appreciated by analysing how different 
de-excitation mechanisms contribute to the residue-production cross section. 



A. Production mechanism in the de-excitation models 

It is instructive to study how the different de-excitation models reconstruct the residue- 
production cross sections as the sum of different production mechanisms. However, we need 
to introduce this discussion by a few important remarks. Firstly, different de-excitation 
models have different, possibly non-overlapping sets of production mechanisms (Sec. 11160 : 
thus, each partition must be seen as model-dependent and cannot be directly compared to 
experimental data or to other partitions. Secondly, although all models internally construct 
some kind of de-excitation-history tree, only a limited, model-dependent amount of infor- 
mation about the decay history is readily available to the user. Figure H] summarises how 
production mechanisms are partitioned in each model; details about each partitioning will 
be given in the model-specific discussions that follow. Each de-excitation mechanism is as- 
signed a colour, which is consistently used in Figs. |5H7l [15] and [181 We attempted to assign 
similar colours to similar mechanisms. In some cases, a specific mechanism in one model 
can be considered equivalent to another mechanism or to the sum of other mechanisms in 
another model. These cases are indicated by arrows and boxes in Fig. HI 

The analysis of the production mechanism will focus on the coupling of the de-excitation 
models with INCL4.5, but the results are qualitatively valid for Isabel, too. 

We start by analysing SMM (Fig. [5]). Firstly, we identify the cascade component of the 
cross section. The rest of the residue-production cross section is partitioned by labelling 
each simulated event with the number nhot of hot fragments that emerge from SMM's initial 
multifragmentation sub-module (and that later de-excite by sequential evaporation). 

Fig. [5] suggests that the plateau cross sections in p+^^^Xe and the IMF cross sections 
in p+^^Fe are almost entirely due to multifragmentation. However, care must be exercised 
with this definition of the multifragmentation contribution. Firstly, the onset of multifrag- 
mentation in SMM is smooth. Close to the multifragmentation threshold, the most probable 
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Figure 4. Summary of the partitioning of production mechanisms in the different de-excitation 
models. Equivalent mechanisms are connected by arrows and are represented by the same or 
similar colours in Figs. [5H71 and [TSVllSI See text for more details about the partitioning. 
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Figure 5. Decomposition of the residue-production cross sections predicted by INCL4.5/SMM accord- 
ing to the number of hot fragments produced in the multifragmentation stage (nhot)- Left panel: 
p+^^Fe. Right panel: p+^^^Xe. Experimental data from Refs. |6|-[9. 



break-up configuration is binary, with one break-up partner much larger than the other. Such 
processes are similar to (and probably indistinguishable from) binary decays, and somehow 
provide a smooth transition to the real multifragmentation regime. Moreover, SMM can pro- 
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Figure 6. Decomposition of the residue-production cross sections predicted by INCL4. 5/ABLA07 ac- 
cording to the de-excitation mechanism. Left panel: p+'^'^Fe. Right paneh p+^^^Xe. Experimental 



data from Refs. 



6-9 



duce events where composite fragments are evaporated during the secondary de-excitation 
of the hot fragments. Whether such events should be counted as multifragmentation is 
unclear. In our analysis, these events are simply classified according to the multiplicity of 
hot multifragmentation products. Thus, the importance of multifragmentation for fragment 
production cannot be easily extracted from the partitioning in Fig. |5]and would require bet- 
ter event labelling, which is unfortunately unavailable at the moment. The predicted cross 
sections for nominal multifragmentation (rihot > 2) , which can be interpreted as upper limits 
for the "real" multifragmentation cross section, are 146.0 mb (^^Fe) and 137.7 mb (^^^Xe), 
which correspond to 18.7% and 10.0% of the respective reaction cross sections. 

Figure [6] displays the partitioning of the INCL4.5/ABLA07 cross sections. Events that 
triggered ABLAOT's multifragmentation module are classified as "multifragmentation"; as in 
the case of SMM, this might include some contamination from secondary fragment evapo- 
ration after nominal multifragmentation. Events that did not trigger multifragmentation 
are catalogued as "fragment evaporation" if one or more fragments were emitted, and as 
"light-particle evaporation" otherwise. Note that events with excitation energies below the 
particle-emission threshold are also classified as "light-particle evaporation". 

One can observe that the cross sections for nominal multifragmentation are much smaller 
than in the case of SMM: 13.3 mb (1.7% of the reaction cross section) for p+^^Fe, and 
4.6 mb (0.3%) for jj+^^^Xe. This is due to ABLAOT's higher multifragmentation threshold. 
According to Eq. ([1]), the freeze-out temperatures for ^^Fe and ^^^Xe are 7.97 and 6.36 MeV, 
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Figure 7. Decomposition of the residue-production cross section predicted by INCL4 . 5/GEMINI++ 
according to the number of asymmetric sphts leading to a given fragment. Left paneh p+^^Fe. 
Right paneh p+^^^Xe. Experimental data from 6|-l9|. 



respectively, which correspond (assuming Ignatyuk's level-density parametrisation [28|) to 
excitation energies of 6.22 and 3.70 AMeV. This should be compared with the typical SMM 
threshold of 3 AMeV. Since the remnant cross section drops fast as the excitation energy 
increases, even a moderate difference in the multifragmentation threshold can result in a 
large cross section difference. 

We anyway stress that most of ABLA07's IMF cross section in p+^^Fe does not origi- 
nate from multifragmentation events. This is at variance with previous claims j^, based 
on INCL4.2/GEM, that the p+^^Fe residue-production cross sections could not be ex- 
plained by evaporation. We stress that, firstly, INCL4.5/GEM provides results similar 
to INCL4.5/ABLA07 on p+^^Fe, and only fails to describe the p+^^^Xe data; secondly, 
INCL4.5/ABLA07 is also able to describe the p+^^^Xe plateau cross sections mostly thanks 
to evaporation alone. Therefore, it is not possible to exclude evaporation solely on the basis 
of the results of one code (GEM) for one system (^^Fe). The details of the evaporation model 
are obviously important, since ABLA07 is able to provide adequate agreement with all the 
experimental data considered so far. 

Finally, the INCL4.5/GEMINI++ partitioning is shown in Fig. [71 De-excitation particles 
are classified according to the number of asymmetric splits that led to their production. For 
particles following from only one asymmetric split, we distinguish if they originated from 
the light or the heavy split partner. For example, if a remnant splits into fragments A and 
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B, with A larger than B, and A subsequently splits into C and D, B will be tallied in the 
"1 asjTiimetric split (light)" histogram, while C and D will be counted as "2+ asymmetric 
splits". Light-particle evaporation {Z < 3) does not influence the count of asymmetric splits. 

Note that we assigned the same colour to the GEMINI++ "0 asymmetric splits" and the 
ABLA07 "light-particle evaporation" components, suggesting that the two mechanisms are 
equivalent (see also Fig. H]). However, emission of Li isotopes in GEMINI++ is described by 
the Hauser-Feshbach evaporation formalism (Sec. Ill B 2|) . and therefore is not counted as 
an asymmetric split. In ABLA07, on the other hand, emission of Li fragments is counted 
as "fragment evaporation". Therefore, the two classes should be considered equivalent for 
emission of fragments with Z > 4. 

The INCL4 . 5/GEMINI++ calculations predict a small symmetric-fission component in 
^_l_i36^g Although ^^^Xe is below the Businaro-Gallone point, a small number of symmetric- 
fission event do occur in remnants with high values of spin and of the Z'^/A ratio. The cross 
section for such events is only about 5 x 10~^ times the reaction cross section and it is not 
visible on the scale of Fig. [71 

Perhaps unsurprisingly, we find that GEMINI++ generates most of the IMF cross section in 
^_l_56pg g^^^^ g^jj ^Yie plateau cross section in p+^'^^Xe through the asymmetric-splitting mech- 
anism. For Z = 3 the contribution with no asymmetric split corresponds to evaporated 
Li nuclei (we remind that emission of nuclei up to Z = 3 is described by the evapora- 
tion formalism; see Sec. Ill B 2p . There is a striking similarity between the right panels 



of Figs. [6] and [71 The contributions of nominal multifragmentation (for ABLA07) and of 
pseudo-multifragmentation 2+-split events (for GEMINI++) in p+^^^Xe are approximately 
of the same magnitude and show a similar Z dependence. This suggests that ABLA07 and 
GEMINI++ also predict similar fragment multiplicities in p+^^^X.e, which is indeed shown to 
be the case in Fig. [HI The curves represent the average number of particles with Z > 3 
as a function of the charge of the fragments that appear in the event. Thus, for example, 
INCL4.5/GEMINI++ predicts that neon fragments {Z = 10) from p^^^^Xe appear in events 
with on average ~ 2.1 particles with Z > 3 (including the Ne fragment itself). Interestingly, 
INCL4.5/GEMINI++ and INCL4.5/SMM predict quite different average fragment multiplicities 
in the ^^^Xe plateau region, despite the residue-production cross sections being very similar. 
In general, the average multiplicity does not seem to strictly correlate with the residue- 
production cross section. This finding manifestly calls for more exclusive observables, such 
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Figure 8. Calculated average multiplicity of particles with Z > 3, as a function of the charge of 
the fragments appearing in the event. See text for more details. Left panel: p+^^Fe. Right panel: 
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as multiplicity distributions and fragment correlations. Such data do exist for p+^^Fe (see 
the following section); the analysis of a SPALADIN-type p+^^^X.e experiment has recently 



been completed and will soon be published 



29| 



As a final remark, we underline that light charged particles emitted during the cascade 
stage might play some role in the determination of cross section for the lowest values of Z. 
By default, INCL4.5 only produces clusters with A < 8, Z < 5. The cascade contribution 
drops off at Z = 4 because none of the selected Z = 5 isotopes has A < 8 (see the beginning 
of Sec. |IV]for the discussion on the FRS isotope selection). Almost 50% of the Z = 3 cross 
section in p+^^^Xe comes from cascade ^'^Li. It is not clear whether heavier clusters might 
significantly contribute to the Z > 4 cross sections, which seem to be slightly underestimated 
by the models. 



V. SPALADIN CORRELATIONS 

We now turn to the analysis of the model predictions for the SPALADIN p+^^Fe data-set 
j^l. Adequate reproduction of inclusive observables is a prerequisite for the study of semi- 
exclusive correlations and/or multiplicity distributions. Therefore, in what follows we will 
only retain the INCL4 . 5 cascade model, which gives residue-production cross sections that 
are in better agreement with the experimental data (see Sec. IIVP . 

The goal of the SPALADIN experiment was to measure observables in coincidence for 
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the 1-GeV p+^^Fe reaction in inverse kinematics. These semi-exclusive measurements were 
obtained at the price of a more sophisticated setup and experimental analysis than typical 
inclusive experiments. In Ref. [o] events generated by cascade/de-excitation models were 
filtered through a Geant4 transport calculation and analysed like the experimental data, 
providing little evidence of convincing multifragmentation signatures. 

Since cascade and de-excitation models are in continuous evolution, it is necessary to pe- 
riodically verify their predictions on sensitive semi-exclusive observables. The INCL4 . 5 and 
GEMINI++ codes, for example, have largely evolved since the publication of Ref. Ig. Moreover, 
the old version of the ABLA code did not allow emission of any nucleus heavier than alpha 
particles, was thus unsuccessful at reproducing inclusive residue-production cross sections 
and was therefore excluded from the study of correlations in Ref. |6|. The new ABLA07 code, 
as proven above, can adequately describe the residue-production cross sections; it is there- 
fore interesting to test its predictions for the SPALADIN correlations. Such considerations 
provide the motivation for this section. 

We start by analysing the decomposition of the residue-production cross sections. Ex- 
perimental events were subdivided into five classes, according to the number of fragments 
{Z > 3) and helium nuclei that were detected. Modelled events were run through a Geant4 
filter that reproduces the experimental setup and simulates detector efficiency, and subse- 
quently categorised just like the experimental events. Figure M displays the result of this 
exercise. Clearly INCL4 . 5/GEMINI++ and INCL4 . 5/ABLA07 yield the best description of the 
experimental data. The appraisal of SMM's results in the IMF region appears to be bi- 
ased by its slight overestimation of the residue-production cross sections. Note also that 
SMM's and ABLA07's overestimation of 3+-fragment events, which are properly reproduced 
by GEMINI++, follow the same trend as the average fragment multiplicities depicted in Fig. [HI 
The INCL4 . 5 / SMM and INCL4 . 5 / GEMINI++ results are quahtatively similar to those presented 
in Ref. O, although both the cascade model and the de-excitation models are sensibly differ- 
ent. 

We now address correlations between fragment charges. We restrict our attention to 
events with at least two detected fragments {Z > 3) and we define Zi and Z2 to be the 
largest and the second-largest observed charges. For these events, we define three bins in 
detected neutron-plus- helium multiplicity (1-2, 3-4 and 5-6), which is expected to be fairly 
correlated with the excitation energy of the cascade remnant. The correlation between de- 
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Figure 9. Partitioning of the 1-GeV p+ Fe cross section according to the number of observed 
helium nuclei and IMFs. Top left panel: experimental data. Other panels: solid lines represent 
the experimental data, histograms represent model predictions. The cascade model is INCL4.5. 
Experimental data from Ref. |6|. 



tected neutron-plus-helium multiplicity and excitation energy was studied with INCL4.2 and 
was shown to be essentially independent of the de-excitation model 30|; however, our results 



contradict this conclusion. Table IIIII shows the average excitation energies per nucleon in 
events with two {Z > 3) fragments as a function of the neutron-plus-helium multiplicity. For 
comparison, we also provide the values computed by Le Gentil et al. with INCL4.2 which 



were claimed to be independent of the de-excitation model [30|. Firstly, we observe that 
the average excitation energies are not independent of the de-excitation model, at least for 
the high-multiplicity bin. Note however that the distributions of excitation energies within 
each bin are broader than the differences in average excitation energies among neighbouring 
bins; thus, the detected neutron-plus-helium multiplicity cannot be interpreted as a precise 
measure of the remnant excitation energy. Secondly, the average excitation energies that 
we find with INCL4.5 are consistently higher than those determined by Le Gentil et al. us- 
ing INCL4.2. Thus, the correlation between neutron-plus-helium multiplicity and excitation 
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provided. See text, Ref. |6| and Ref. |30| for more details. 

Table III. Remnant excitation energies per nucleon in 1-GeV p+^^Fe in events with two detected 
fragments with Z > 3, for three bins in detected neutron-plus-helium multiplicity, as calculated 
by the three de-excitation models used in this work (coupled to INCL4 . 5) and by the INCL4 . 2 
calculations performed by Le Gentil et al. 6|. The values should be interpreted as mean value it 
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Figure 10. Distributions of Zi — Z2 (see text) for three different bins in neutron-plus-helium 
multiplicity, in events with at least two fragments with Z > 3. See the text for more details about 
the excitation energies sampled in these bins. The cascade model is INCL4.5. Experimental data 



from Ref. 



energy depends at least on the cascade model; neutron-plus-helium multiplicities cannot 
represent a universal, model-independent measure of the remnant excitation energy. 

Figure [10] shows the distributions of the Zi — Z2 difference in the three multiplicity 
bins. As explained in Ref. [g], the experimental data indicate that asymmetric {Zi ^ Z2) 
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Figure 11. Average remnant excitation energy as a function of the value of Abound, which is defined 
as the sum of all the detected charges greater than one. The error bars represent the width of the 
excitation-energy distributions for INCL4 . 5/GEMINI++. The black line was extracted from Ref. 
and was computed with INCL4 . 2/GEMINI. 



and symmetric {Zi ^ Z2) charge configurations are favoured at low and high excitation 
energy, respectively. All the models reflect this qualitative trend, although only GEMINI++ 
and ABLA07 are able to reproduce the absolute cross sections with good accuracy. SMM's 
overestimation is qualitatively consistent with its predictions of IMF yields; in fact, we 
observe that the cross sections calculated with a given model, if summed over Zi — Z^ 
and over the multiplicity bins, are quantitatively comparable to the IMF-production cross 
sections shown in Fig. O 

Compared to Le Gentil et al.'s results, our work confirms that INCL4.5/GEMINI++ pro- 
vides the best description of the Z\ — Z2 distributions. The new ABLA07 model, when coupled 
with INCL4 . 5, also provides a very good reproduction of the experimental data. However, 
one should also remark that the shapes of the INCL4 . 5/SMM distributions are now quite sim- 
ilar to the experimental data, which was not the case in Ref. [g, and only the normalisation 
seems to be consistently off by a factor of about two. This is quantitatively consistent with 
INCL4.5/SMM's overestimation of the residue-production cross sections (Fig. |2]) and raises 
an interesting question, i.e. whether it might be possible to adjust INCL4.5/SMM to better 
reproduce residue-production cross sections and the SPALADIN observables at the same 
time. 
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We finally turn to fragment multiplicities. We first need to define the Abound variable as 
the sum of all the detected charges with Z > 2. This quantity was previously found to be 
negatively correlated with the excitation energy of the cascade remnant, and the correlation 
was found to be independent of the de-excitation model Q|. Fig. [TT] demonstrates that 
both these properties stay true when the calculations are performed with INCL4 . 5. The E*- 
■Z^bound correlation, however, is found to be slightly different for INCL4.2 and INCL4.5. This 
indicates that Abound should not be considered as a universal, cascade-model-independent 
measure of the excitation energy, although it seems relatively robust with respect to the 
choice of the de-excitation model. 

The average fragment multiplicities (i.e. the average number of fragments with a given 
charge produced in a reaction) are plotted in Fig. [121 fo^' fragment charges between 2 and 
7, as functions of Abound- All multiplicities rise as 2'bound decreases, reflecting the positive 
correlation between fragment multiplicity and excitation energy. Note that production cross 
sections for Z = 2 were measured by the SPALADIN collaboration, but they were not 
plotted on Fig. [2] because we summed the calculated isotopic cross sections over the isotopes 
measured in the FRS experiments (see Sec. II Vp . 

The first striking result is that, for a fixed fragment, all the curves have a very similar 
shape. If Abound is interpreted as the excitation energy of the cascade remnant, this indicates 
that all the de-excitation models predict a similar dependence of the fragment-emission 
probability on the excitation energy. The overall level of the curve is charge-by-charge 
correlated with the IMF yields in Fig. [2l Thus, for example, ABLA07 and SMM predict 
too large multiplicities and too large production cross sections for Z = 3. Besides this 
correlation, it is not clear whether there is a lesson to be learnt from these observables. The 
strong even-odd staggering at low Abound is reproduced by all models. 

Summarising, the quest for model-independent measures of the remnant excitation energy 
is still open. We checked that neutron-plus-helium multiplicities and the Abound variable 
are at least sensitive to the cascade model. The new GEMINI++ and SMM versions yield 

n 

predictions similar to those reported by Ref. [g, although the cascade model used in the 
present work yields rather different excitation-energy and remnant-mass distributions. The 
ABLA07 and GEMINI++ models can reproduce most of the considered observables. Given 
the very small fraction of nominal multifragmentation events predicted by INCL4.5/ABLA07 
(Sec. IIV A]) , we confirm that explanation of the SPALADIN data does not require any strong 
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Figure 12. Average fragment multiplicities as a function of the Abound variable (see text). The 
cascade model is INCL4.5. Experimental data from Ref. [6. 



multifragmentation component . 



To our knowledge, no published correlation data exist for the p+^^^Xe reaction around 1 
GeV. However, a SPALADIN-type experiment was performed in April 2009 and the results 



of the analysis are due to be published soon 



29|. 
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VI. LONGITUDINAL- VELOCITY DISTRIBUTIONS 



Insight about the de-excitation mechanism can also be gained by examining the kine- 
matics of the decay products. Sequential binary splits should produce kinematical patterns 
reminiscent of the decay barriers; multifragmentation, on the other hand, is expected to pro- 
duce fragments with broad, structureless velocity distributions. We discuss here the velo city 
distributions measured in the context of the FRS experiments considered in Sec. IIVI j^, llo| . 
The emission velocities of fragments from the 1-GeV p+^^Fe and p+^^^Xe reactions are mea- 
sured using forward spectrometry techniques. Reactions are studied in inverse kinematics, 
i.e. as a 1-AGeV ^^Fe or ^'^^Xe beam impinging on a target. Most of the de-excitation 
products are focalised in a cone around the beam axis. The experimental angular acceptance 
somewhat depends on the trajectory azimuth, but it is on average equal to 15 mrad. Parti- 
cles that satisfy the acceptance cut are identified by mass and charge and their longitudinal 
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velocity (the component of the velocity along the beam axis) is measured. Refs. |9| and 
report measured longitudinal- velocity distributions (LVDs) for several nuclides with A> Q. 

We used our cascade/de-excitation tools to calculate the longitudinal- velocity distri- 
butions for the same nuclides. A detailed comparison with the experimental LVDs re- 
quires knowledge of the azimuthal dependence of the FRS angular acceptance and a three- 
dimensional macroscopic transport calculation of the reaction products in the spectrometer. 
We limited ourselves to a simpler approach: we assumed the angular acceptance to be inde- 
pendent of the trajectory azimuth and equal to the experimental average value of 15 mrad. 
In other words, the acceptance of our simulation is a circular cone in velocity space, centred 
on the beam axis, with vertex in the origin and aperture of 15 mrad. This prevents a refined 
quantitative comparison of our results with the experimental data, but the emerging trend 
is nonetheless clear, as we will show in the following. 

Figure [13] shows the calculated LVDs for ^^Fe+^H. Each distribution is separately nor- 
malised to one. Note that, in this and all the following figures, the longitudinal velocities 
refer to the rest frame of the ^^Fe projectile, with the proton impinging with negative velocity. 
This choice was made for consistency with the experimental data plotted in Fig. [TH 

We note that a few short-lived nuclides are present in the SMM results (^Li, ^Be and ^B). 
These nuclides would typically decay before being detected by the experimental apparatus. 
The decays could in principle populate other IMF species and modify their LVDs, but the 
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Figure 13. Calculated longitudinal-velocity distributions for 1-AGeV ^^Fe+^H in the beam rest 
frame. Red solid line: INCL4 . 5/GEMINI++; green dashed line: INCL4. 5/ABLA07; blue dotted line: 
INCL4.5/SMM. 



nuclides above entirely decay in nucleons and alpha particles. Therefore, we can neglect 
them in the following discussion. 

One observes that all models produce similar, single-peaked distributions for the heaviest 
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Figure 14. Experimental longitudinal-velocity distributions for 1-AGeV ^^Fe+^H in the beam rest 
frame. Adapted from Ref. 0. 



IMFs (say for A > 9). Only for the Ughtest IMFs can we observe differences among the 
model predictions, with GEMINI++ and ABLA07 often producing double-peaked distributions, 
whereas SMM typically yields fiat distributions. These predictions should be compared with 
the measured distributions jo]. Fig. 10], which are reported in Fig. [14] for convenience of 
the reader. In stark contrast with Napolitani et al.'s claims, we find that binary decay 
does not imply sharp Coulomb holes in the velocity distributions. Indeed, the shapes of 
the measurements distributions seem to be best described by GEMINI++; compare e.g. the 
double-peaked structure of the ^'^Li distributions, where the Coulomb peaks predicted by 
GEMINI++ are possibly even too weak to account for the measured shape. 

How can GEMINI++ produce single-peaked LVDs by relying on its binary-decay mecha- 
nism? This question can be answered by partitioning LVDs according to the production 
mechanism, as done in Fig. [T^] Firstly, we observe that nuclei that follow from two or 
more asymmetric splits (cyan component) expectedly produce single-peaked distributions. 
The "0 asymmetric splits" component, which is only present in Li isotopes, corresponds 
to direct evaporation and shows a sharp Coulomb hole. However, contrary to what one 
would naively expect, nuclei following from one asymmetric split produce only mildly struc- 



26 




Longitudinal velocity [cm/ns] 

Figure 15. Decomposition of simulated longitudinal-velocity distributions in ^^Fe+^H for some 
selected nuclides, according to the partitions summarised in Fig. 



tured LVDs. Distributions associated with light split partners do show Coulomb holes, 
although secondary de-excitation (evaporation of light particles) after the split somewhat 
blurs the peaks. However, and more importantly, the heavy partner of asymmetric splits 
typically picks up very little recoil, but can retain enough excitation energy to lose much 
of its mass and be eventually detected as an IMF with small longitudinal velocity. Thus, a 
sizable fraction of the single-peaked contribution to the GEMINI++ LVDs in Fig. [15] comes 
from the "1 asymmetric split (heavy)" (dark-blue) mechanisms, which we may term the the 
de- excitation-residue component. This result somewhat contrasts with the interpretation 
suggested by SMM, which entirely attributes the single-peaked component to nominal multi- 
fragmentation events. Finally, observe that the ABLA07 "fragment evaporation" component 
(blue), which contains both evaporated fragments and fragment-evaporation residues, shows 
no stark Coulomb structure for A > 10, consistently with the conclusions suggested by 
GEMINI++. 

We stress that the participation of de-excitation residues to IMF distributions is possible 
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only because cascade remnants in p+^^Fe are relatively close in A and Z to the IMF mass 
region (Fig. [1]). The residue component somewhat provides a background noise over which 
true multifragmentation signatures are superimposed. One of the motivations for studying 
the p+^^^Xe reaction is exactly that the residue component is expected to be completely 
negligible in the IMF region. It is therefore of great interest to consider how well GEMINI++ 
performs at reproducing the measured LVDs for this reaction. 

As displayed in Figs. [16] and [T71 the LVDs predicted by GEMINI++ show clear Coulomb 
holes over all the considered mass and charge range. Indeed, the residue component is 
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negligible and GEMINI++ is not able to reproduce the experimental distributions 
which become rather fiat from Z ~ 6, as in the case of ^^Fe. Surprisingly, however, even 
ABLA07 and SMM predict double-peaked LVDs that are very similar to GEMINI++'s, even for 
relatively heavy nuclides such as the Si isotopes (Fig. [T71) . No model seems to be able to 
account for the single-peaked component that clearly dominates the experimental LVDs for 

This unanticipated results clashes with the widespread belief that multifragmentation 
should yield single-peaked velocity distributions, especially for high fragment multiplici- 
ties. SMM accounts for the exact propagation of the hot multifragmentation products in 
their mutual Coulomb field, and is therefore supposed to yield reasonable predictions of the 
asymptotic velocities, if the initial conditions are realistic. Fig. [TS] demonstrates that this 
mechanism does not yield single-peaked LVDs. Even many-body (> 3) break-up configura- 
tions (dark-red component) bear clear signs of Coulomb repulsion. This suggests that most 
of the break-up configurations must be quasi-binary, with one or two large fragments com- 
pletely dominating the Coulomb dynamics, possibly accompanied by nucleons and very light 
charged particles. Note that ABLAOT's multifragmentation mechanism seems to produce flat- 
ter distributions, although its cross section is largely insufficient to explain the experimental 
shapes. 

Finally, we comment briefly about the contribution of dynamical emission of cascade 
clusters to the LVDs. We remind the reader that INCL4 . 5 by default produces clusters 



12j . Indeed, Figs. [T5] and [TSl show a 



through a coalescence mechanism up to A = 8 included 
cascade contribution for the lightest IMFs. Cascade clusters, which are high-energy particles 
in direct kinematics, constitute an asymmetric tail that extends in the backward direction 
in inverse kinematics. The cascade-cluster tail is sometimes responsible for a large forward- 
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Figure 16. Same as Fig. m for 1-AGeY ^^exe+iH. 



3ackward asymmetry of the LVD. A similar signature was observed in the ^^^Xe+^H data-set 



10|, but was attributed to fluctuations in the recoil momentum of the cascade remnant. The 



INCL4.5 model suggests a different interpretation. 
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Figure 17. Same as Fig. [T6l for 9 < Z < 14. 

VII. TIME INTERVAL BETWEEN FRAGMENT EMISSIONS 

We have so far considered residue-production cross sections, correlations among de- 
excitation products and distributions of longitudinal velocity of the emitted fragments. For 
most of these endpoints, the INCL4.5/GEMINI++ model provided the most accurate descrip- 
tion of the experimental data (with the exception of the p+^^^Xe LVDs, which no model 
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Figure 18. Same as Fig. [El for "^Xe+^H 



seems to be able to reproduce). This result can be taken as evidence that a multifragmen- 
tation model is unnecessary to describe the reactions studied in the present paper. The 
solidity of this argument, however, relies on the internal consistency of the application of 
the INCL4.5/GEMINI++ model to the systems considered. We will now proceed to show that 
the time interval between fragment emissions for highly excited cascade remnants becomes 
comparable to the typical multifragmentation timescale. 

The GEMINI++ model keeps an internal clock of the decay process, which is readily avail- 
able to the user. For a compound nucleus with decay width F, the decay time is sampled 
from an exponential distribution with time constant U/T. Thus, we select events with two 
or more asymmetric splits {two-split events) and compute the interval length Atgpiit between 
the earliest and the second-earliest emission. Note that this definition is the event-based 
equivalent of the "2+ asymmetric splits" particle classification above (Sec. IIV Ap . If more 
than two splits occur during an event, we only consider the two earliest; the idea is that 
we are interested in short intervals and emissions become more separated in time as the 
excitation energy is evacuated. Note that not every asymmetric split leads to observed 
fragments; if the excitation energy of the emitted fragment is sufficient, it can completely 
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Figure 19. Distributions of time intervals between asymmetric splits. Curve colours, from blue to 
red, correspond to increasing values of the excitation energy. Each curve is normalised to one. Left 
panel: p+^^Fe. Right panel: p+^^^Xe. 

disassemble in light charged particles and remain unobserved. Thus, two asymmetric splits 
do not necessarily correspond to three observed fragments. 

We then define twelve bins in excitation energy per nucleon by requiring that two-split 
events be uniformly partitioned over the bins. Since there are only few two-split events 
at low excitation energy, the first bin is very broad (from zero to 4.2 AMeV for ^^Fe and 
2.8 AMeV for ^^^Xe). Finally, for each bin we construct a distribution of interval lengths. 
This procedure permits studying how time intervals evolve as a function of the excitation 
energy of the cascade remnant. 

Fig. [19] displays the distributions of time intervals for the twelve excitation-energy bins. 
The lowest excitation-energy bins are given above. The highest excitation-energy bins range 
from 8.8 AMeV (^^Fe) and 4.6 AMeV (^^^Xe) to infinity. One notices that the distributions 
corresponding to low excitation energies are close to exponential. All distributions are 
very broad, with heavy tails extending up to several thousands fm/c; for this reason, it 
is inappropriate to characterise the distributions using their mean and/or their standard 
deviation. In what follows, we will rather rely on the median and the interquartile range, 
which are more robust. 

The dotted lines in Fig. EH] represent the energy-differential cross section for two-split 
events (right scale). One immediately notices that two-split events are concentrated at 
higher excitation energies per nucleon in ^^Fe than in ^^^Xe. The ^^Fe distribution extends 
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Figure 20. Median (solid lines) and quartiles (dashed lines) of the distributions of time-interval 
lengths between fragment emissions (Aigpiit), as a function of the initial excitation energy per 
nucleon. The model is INCL4. 5/GEMINI++. Red (blue) lines correspond to the p+^^^Xe (p+^^Fe) 
reaction. Cross sections for two-split events are superimposed as dotted lines (right scale). 

up to very high excitation energies, comparable to or larger than the total binding energy of 
the remnant. Such remnants are rare but nevertheless possible. Interestingly, the maxima 
of the two-split cross sections are located at about 5.8 (^^Fe) and 3.6 AMeV (^^^Xe), which 
are similar to ABLA07's multifragmentation thresholds cited in Sec. IIV A[ This strengthens 
the (perhaps coincidental) similarity between ABLAOT's multifragmentation and GEMINI++'s 
two-split events that was observed in Sec. IIV A[ 

On the left scale of Fig. |20]we report the medians (solid lines) and quartiles (dashed) of the 
interval-length distributions in each excitation-energy bin. The fact that the quartile curves 
are approximately parallel to the median curve indicates that, up to a scale factor (note the 
logarithmic scale on the Atgpiit axis), the shapes of the distributions are approximately the 
same in all the bins. If we take the median as a measure of location, we can observe that 
typical interval lengths decrease as the initial excitation energy increases. Note also that, for 
the same excitation energy per nucleon, interval lengths are larger for ^^Fe than for ^'^^Xe. 

INCL4 . 5/GEMINI++ thus predicts that typical intervals between fragment emissions last 
a few hundreds fm/c at the onset and reach 70-75 fm/c at the peak of the two-split cross 
section (independently of the target). For higher excitation energies, even shorter times are 
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expected. These numbers are comparable to the typical multifragmentation timescale of 



a few tens of fm/c 
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32|, suggesting a continuous transition between sequential binary 



decay at low energy and the expected multifragmentation regime at high energy. This 
aspect suggests that a binary-decay model can generate final states that similar to those 
that a multifragmentation model would produce. 

The time interval between fragment emissions represents an upper bound for the interval 
between any two consecutive binary decays. As the excitation energy increases, this time 
eventually becomes comparable to the relaxation time of the system. Under these condi- 
tions it is difficult to justify the compound- nucleus hypothesis, which assumes a completely 
equilibrated system. However, it is difficult to provide quantitatively accurate estimates of 
the relaxation time of a highly excited nuclear system. Besides equilibration, however, the 
asymptotic (observable) escape velocities of the emitted charged particles are sensitive to the 
length of the interval between emissions 31| and should in principle be determined from the 
solution of the equations of motion of the emitted fragments in their mutual Coulomb field. 
In GEMINI++, as in most statistical de-excitation codes, it is assumed that decay products 
have already attained their asymptotic velocity before they undergo any subsequent decay. 
The importance of an exact solution could then be evaluated by studying observables that 
are sensitive to the de-excitation kinematics, such as LVDs. Note however that SMM does 
include a numerical solver for the Coulomb trajectories of the hot fragments, but it is still 
unable to reproduce the experimental LVDs for p+^^^\e. 



VIII. CONCLUSIONS 

We have used the tools of coupled intranuclear-cascade and statistical-de-excitation mod- 
els to search the 1-GeV p+^^Fe and p+^^^Xe reactions for signatures of multifragmentation. 
The choice of the cascade model has some influence on the distributions of remnant size and 
excitation energy; in particular, dynamical emission of clusters during the cascade stage has 
a sensible influence on the remnant-mass distribution and, thus, on the residue-production 
cross sections. This leads to a rather large sensitivity of calculated residue-production cross 
sections on the cascade model. For the purpose of this study, we chose to flx the cascade 
model by requiring that it correctly reproduce residue-production cross sections close to the 
target nuclide, which are typically understood as due to the evaporation of lowly-excited 
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cascade remnants and are rather insensitive to the choice of the de-excitation model. 

Calculations indicate that the inclusion of a multifragmentation stage is not crucial for 
adequate prediction of residue-production cross sections. We thus confirm the insensitiveness 
of this observable to the de-excitation mechanism. However, different de-excitation models 
propose widely different reconstructions of the residue-production cross sections in terms 
of elementary processes, suggesting that semi-exclusive observables can help discriminate 
among different de-excitation mechanisms. Comparisons with measured fragment-helium 
correlations, Zx — Z2 distributions and IMF-gated Abound distributions favour binary 
de-excitation models such as GEMINI++, or models predicting very small multifragmentation 
cross sections such as ABLA07. We conclude that a multifragmentation model is not necessary 
for the description of inclusive and semi-exclusive observables. This does not mean that the 
presence of multifragmentation is ruled out in these reactions, but rather that binary decay 
can generate final states similar to those produced by multifragmentation models, at least 
close to the multifragmentation threshold. 

Somewhat ambiguous conclusions can be drawn from the qualitative study of longitudinal- 
velocity distributions. Contrary to previous claims jol, we find that pure binary decay can 
account for the distributions measured in p+^^Fe. The observed single-peaked component 
can be ascribed to multifragmentation, de-excitation residues, or both, depending on the 
de-excitation model considered. On the other hand, none of the considered de-excitation 
models can explain the existence of the observed single-peaked component in jj+'^^^Xe. Even 
numerical integration of the Coulomb trajectories of the multifragmentation products, as 
implemented in SMM, predicts double-peaked longitudinal-velocity distributions. Therefore, 
the shape of the longitudinal- velocity distribution is an ambiguous signature of de-excitation 
mechanism. 

The INCL4.5/GEMINI++ calculations suggest that ^j+^^^Xe residues with 4 < Z < 40 
are mostly produced in events with one asymmetric split, with no contribution from de- 
excitation residues; it is far from obvious that a similar mechanism can produce single- 
peaked longitudinal-velocity distributions. If the single-peaked component in p+^'^^Xe must 
be ascribed to multifragmentation, we would expect the multifragmentation signature to be 
even more visible in p+^^Fe; but that does not appear to be the case. 

We have also studied the time interval between asymmetric splits in INCL4.5/GEMINI++. 
At the highest excitation energies per nucleon, the model predicts interval lengths compa- 
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rable with the typical multifragmentation timescale. This again suggests a smooth tran- 
sition between the binary-decay and the muhifragmentation regimes and illustrates how 
binary decay can generate multifragmentation-like final states, as mentioned above. It is 
not clear whether equilibration times of the order of the time interval between asymmetric 
splits are sufficiently long to justify GEMINI++'s compound-nucleus hypothesis; in any case, 
closely-packed binary emissions of charged fragments are expected to distort the asymptotic 
Coulomb velocities. This effect is not accounted for in any of the models considered in the 
present work. 

In conclusion, binary decay yields a satisfactory description of most of the observables 
considered in this paper. The application of binary-decay models to cascade remnants with 
very large excitation energies generates final states that resemble those produced by mul- 
tifragmentation models. The good agreement of INCL4.5/GEMINI++ and INCL4.5/ABLA07 
with the experimental data considered in this paper probably indicates that events with very 
high excitation energy per nucleon do not significantly contribute to the studied observables. 
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